UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



LIBRABY ROUTINE H 3 - SO 



TITLE 

TYPE 

NUMBER OF WORDS 

TEMPORARY STORAGE 



ACCURACY 
DURATION 



PRESET PARAMETERS' 



ERTRY 



Minimisation of a function of Tvo Variable (DO! or SADOI ) 
Closed 

3 words at S3 

3 words at S^ 
3 words at Sp 

Depends upon the condition of the function. 
A "minor" cycle takes T = (Ut) ms where n Is the number 
of "minor" cycles in the "major cycle". Add 1 ms for each 
decrease in mesh size. 

3 00 F 00 aF where a, a+1, a+2 are the addresses of 

f(xy), x, y on entry and exit, 

b, b+1, b+2 are the addresses of 5, 6 , 5 , 

c, c+1, c+2 are working spaces , 
a is the factor by which the mesh size is 
decreased, 

e is an end constant such that only mesh 
sizes larger than 6 will be used, 
r is the address to the left hand side of 
which control is transferred after the 
best values using a given mesh size have 
been obtained. 

s is the address of a closed subroutine 
which takes x from c+1 (1S5), J from 

c+2 (2S5) and places f (xy) in R.^ 
Enter with initial values of x, y, 5 in 1S3, 2S3, S^, 
respectively, and 

5Q pF 



k 


00 T 00 bF 


5 


00 F 00 cF 


6 


00 F 00 OJ 


7 


00 F 00 eJ 


8 


00 F 00 rF 



9 00 F 00 sF 



p+1 



q = address of this program 



26 qF 

Routine will be left with best values of f (xy), x, y 
in S3, 1S3, 2S3, respectively. 
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DESCRIPPION Define 



f x = f (x - S,y) - f(x + 5,y) 
f = f(x,y - 5) - f(x,y + 5) 

1 - The routine forms f and f 

x y. 



2 - 5 = 
x 



1/2 f 5 
' x 



1/2 |f x | + 1/2 |f y | + 2~ 59 



1/2 f 5 
6 = _ l_i_ 



y. 



1/2 |f x | + 1/2 |f y | + 2* 59 



are formed. Hates The initial mesh size should be chosen 

so that f and f do not exceed capacity. 
x y . J 

3 - f (x + 5 , y + 5 ) is formed,, and tested against the 

previous f(x,y). If the function has decreased f (x,y) 

is replaced by f (x + 5 , y + 6 ). The process is then 

x y 

repeated with the same 5 and 5 « This constitutes a 
"minor cycle". 

k - When the above process ceases to improve the function, 

we return to steps 1 and 2 and compute a new 6 and 5 « 

x y 

Steps 1-4 constitute a "major cycle". 

5 - When the first " min or" cycle of a given "major" cycle 
fails to improve the function,, indicating that no further 
improvement is to be expected with the current mesh size, 
control is transferred to the left hand side of the address 
in S8. The programmer may here place any program desired ': 
e.g. printing of intermediate results (best curreat f (x,y), 
x, and y are always in S3* 1S3* 2S3)j alteration of the 
program depending upon possible knowledge of the function^ 
etc. The program inserted in S8 should terminate by trans- 
ferring control to the left hand side of (q + hi)* If 
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the programmer does not wish to leave this code, 00 F 
00 .(q+UljF should be entered in S8 during read- in. 

6 - Upon returning to the code, 5, the mesh size will be 
replaced by a 6 where a is entered in S6 during read-in. 
The steps 1-5 are then repeated. 

7 - The steps 1-6 will be repeated (N+l) times where N is 

N 
the largest number such that a 5 < e. € is to be entered 

in ST during, read -in. 

NOTES 1 - x , y Q , b The initial, guesses for s, y, 8 are to be 

stored in 1S3, 2S3, and Sk "before entering this subroutine. 

2 - A function, poorly conditioned in having a much larger 
gradient in one direction (s&y x) than in the other (say 
y), can deceive this routine* Such deception can sometimes 
be penetrated by re-entering the code with a larger 

mesh size after the code ha8 Jaeen left, retaining the 
best x and y (in 1S3 and 2S3| Obtained on the previous 
entry. 

3 - By properly programming tw£ subroutine to compute 

f (xy) to neglect y, the routi^ may be used to minimize 
a function of one variable. 
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DAT E July 20^ 195^ rt.j2/l6/58 
CODED B Y J. H.Snyder 
APPROVED B Y Q^^^J^O 



±gr 



- 
LOCATION 


ORDER 


NOTES PAGE 1 


.0 


00 K(H3) 

S5 F 

L4 9L 




Plant link 


1 


42 43L 








L5 1S3 




Xo y to 1S5, 2S5 


2 


40 1S5 
L5 2S3 




; 


i 

3 


40 2S5 
50 3L 




! 


4 


26 S9 
40 S3 




j 


5 


L5 HL 


From 35 


Set switch for failure on first try 




42 35L 


and 43 




6 


L5 1S3 
L4 S4 






7 


40 1S5 
L5 2S3 






8 


40 2S5 
50 8L 




i 
I 


9 


26 S9 


• 


1 

i 

i 




10 IP 




l/2 f (x + 5, y) to 1S4 \ 


10 


40 1S4 
L5 1S3 






11 


LO S4 
50 S8 






12 


40 1S5 
50 12L 




! 
1 

i 


13 


26 S9 




1 




10 IF 




l/2 f(x - §, y) in r" j 


14 


LO 1S4 




. 




40 1S4 




1/2 f in 1S4 

' X 


15 


L5 1S3 
40 1S5 




• 

i 


16 


L5 2S3 




1 

i 




L4 S4 





LOCATION 


ORDER 


Ntffes PAGE 2 | 


H 


IT 


1*0 2S5 
50 17L 




■ 

l/2 f to 2Sl* 




18 


26 S9 
10 IF 




1 




19 


40 2Sl* 
L5 2S3 




i 

i 

1 




20 


LO Si* 
50 5L 








21 


40 2S5 
50 21L 




! 

i 




22 


26 S9 
10 U 




i 

s 




23 


LO 2Sl* 
1*0 2Sl* 


• 


1 
s 




2k 


19 38F 
L6 1S1* 




i 




25 


L6 2Sl* 
1*0 S5 




' 




26 


L5 ISl* 
66 S5 




. 




27 


7J Si* 
1*0 isi* 




5 in lSl* 

X 




28 


L5 2Sl* 
66 S5 








29 


7J Si* 
1*0 2Sl* 




5 in 2Sl* 

X 




30 


L5 1S3 
Ll* lSl* 


From 1*0 






31 


1*0 1S5 
L5 2S3 




y + 5 in 2S5 




32 


Ll* 2Sl* 
1*0 2S5 




y+5 v ia 2S5 




1 35 


50 33L 




Waste 




I 


50 35L 


j 







LOCATIOH 



ORDER 



34 



35 



56 



37 



38 



39 



40 



41 



42 



43 



44 



45 



26 S9 
40 S5 
LO S3 
36 ( )F 
L5 S5 
40 S3 
L5 1S5 
40 1S3 
L5 2S5 
40 2S3 
L5 20L 
42 35L 
26 30L 
00 F 
50 44L 
TJ S4 
40 S4 
LO 45L 
36 5L 
22 ( )F 
00 F 
00 S6 
00 F 
00 ST 



By 5 and 39 



By 2 






f (x + 8 . y + 



) in R. and S5 
■' 1 



f - f (before ) in R 1 

If negative, success 
I Replace f Q , x Q , y Q and set 

switch for failtire after first try 



and 
repeat 

Waste 



8 i + i * h 

and S^ 
6 i+l " £ ta E l 



Link 
5 



i+1 
6, 



= a 



= € 



